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We apply a recently developed phenomenological theory of polar liquids to calculate the repulsive 
pressure between two hydrophilic membranes at nm-distances. We find that the repulsion does 
show up in the model and the solution to the problem fits the published experimental data well 
qualitatively and even quantitatively. Moreover, we find that the repulsion is practically independent 
of the temperature, and thus put some extra weight in favor of the so called hydration over entropic 
hypothesis for the membranes interactions explanation. The calculation is a good "proof of concept" 
example a sufficiently sophisticated continuous water model application to a non-trivial interaction 
on nm-size objects in water arising from long-range correlations between the water molecules. 



The solvent induced forces play an important role in 
Nature, key nano- and bio-technological applications, 
drive various self-assembly phenomena in cells and mem- 
branes p], vesicle- membrane, and -synapse fusion [21 E], 
protein folding [4]. Understanding of these phenomena is 
crucial in a number of practically important applications 
such as drug design, drug transport optimization, and 
the design of nano-particles drug transport and delivery 
systems [S]- 

One of the most conceptually simple while still an im- 
portant and an easily experimentally observable example 
of such interactions is provided by the short range re- 
pulsion between phospholipid membranes, originally dis- 
cussed in |6] [3- Since the original discovery there has 
been a lot of efforts to understand the nature of the re- 
pulsive forces, as described in [J and the refs. therein. 
Most of the time the interactions are analyzed within 
some kind of a two-body approximation in such a way 
that the total water mediated pressure between the two 
parallel lipid membranes is represented as a the sum of 
the two components: the direct pressure, arising from the 
direct interactions between the opposing membranes, and 
the hydration component, associated with the interaction 
of the membranes with the intervening water molecules 

Since the membranes are hydrophilic, the exclusion of 
the water leads to the energy loss, and hence on the 
molecular level the hydration pressure can be associ- 
ated with the orientation polarization of water near the 
membrane surfaces |9][10l|TlJ. An alternative approach 
summarized in [12] suggests the dominant role the water 
molecules ordering next to the surfaces and thus empha- 
sizes the water entropy contribution to the free energy of 
the system. None of the effects are easy to grasp within 
any kind of a simple continuous water model, which are 
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normally designed to reproduce the effects of continuous 
electrostatics in polar solutes, e.g. |131 [Ti] 115] . The rea- 
son is that due to long-range electrostatic interactions 
between the molecular dipoles the correlations in water 
are collective and survive at very long distances up to 
Inm. Therefore the effects of the ordering may lead 
to cluster formations and phase transitions phenomena 
[m [m UHl [m EQ] EH Eg, appearance of strong and 
sufficiently long range interactions of non-electrostatics 
nature [UESjEl]- 

The alternative to the continuous solvation models is 
Molecular Dynamics (MD) |25j of the body of inter- 
est immersed in a tank of water molecules in a realis- 
tic force field or even within quantum mechanical set- 
tings. Though such an approach may in principle pro- 
vide ultimately accurate predictions, the calculations are 
time consuming and pose a number of challenges stem- 
ming, e.g. from long relaxation times of water clusters. 
One possible way to bridge the "simulation gap" is to 
develop advanced of continuous solvation models, such 
as |22l [23l [26] [27] and test its limits to make sure the 
models include realistic and important interactions. In 
what follows we take the polar liquid phenomenology and 
calculate the repulsive pressure between two hydrophilic 
membranes at nm-distances. We find that the repul- 
sion does show up in the continuous model and that the 
solution to the problem fits the available the published 
experimental data well both qualitatively and quantita- 
tively. The pressure turns out to be practically indepen- 
dent of the temperature, which strongly supports the hy- 
dration nature of the water-assisted interactions between 
the membranes. The interaction appears to be generic, 
should arise between any hydrophilic bodies and vanishes 
as soon as the objects in question are separated farther 
then Rd ~ Inm, the characteristic polarized water do- 
main size identified in the model. 

The continuous polar liquid model introduced in [23] 
is capable of describing both the short- and the long- 
range features of a polar liquid in a single theoreti- 
cal framework. Originally it was applied to calculate 



water-assisted interactions of macroscopic bodies with 
hydrophobic interfaces of various shapes and charges. It 
was shown that the competition between the short range 
(hydrogen bonding) and the long-range dipole-dipole in- 
teractions of the solvent molecules leads to appearance of 
strong, long range and orientationally dependent interac- 
tions between the objects, which, in principle, can be re- 
sponsible for various self-assembly processes in biological 
systems. Within the suggested model a polar liquid it- 
self is characterized by a complicated fluctuating thermal 
state, ordered at sufficiently short scales within a single 
domain, and completely disordered at larger distances. 
This physical picture has far reaching consequences, es- 
pecially at solvent-solute surfaces, where the ferroelectric 
film of solvent molecules may be formed [28]. In accor- 
dance with the MD simulations |20l |2T] the vector model 
predicts the BKT-like topological phase transitions at 
solvent-solute interfaces |291 130] . Most of our earlier re- 
search was confined to interactions of hydrophobic bod- 
ies. In this Paper we try the model and provide the solu- 
tion to the phospholipid membranes repulsion problem, 
which is a hydrophilic bodies interaction example. 

The model description of a polar liquid, e.g. water, 
proceeds as follows. Each of the molecules within the liq- 
uid {j = 1,2, . . . , N) is assumed have a vector property: 
the static electric dipole moment — doSj, where dp is 
its magnitude, and Sj is the unit vector in expressing the 
orientation of the molecule. Having this in mind it is pos- 
sible to develop a vector field theory in which the liquid 
is described by a local mean value of the molecular po- 
larization vector s(r) — {d)/dQ, where (d) = (d^) stands 
for the statistical average of molecular dipole momenta 
d over a small but sufficiently large volume of the liquid 
containing macroscopic number of molecules. Accord- 
ingly, the model expression for the free energy functional 
is given by: 

a.,P—x,y,z 

+ I dVV{s^) + 1 dV^El-J dV^P(r)Ee(r). (1) 

The polarization vector of the liquid, P = PqS where 
Pq — nodo, no is the molecular density, and Eg = — V(/?e 
is the external electric field in the absence of the liquid. 
Similarly, Ep — —Wipp is the polarization electric field, 
produced by the polarization charges within the liquid 
characterized by the polarization charge density pp = 
— divP. The polarization potential (pp should be found 
from the Poisson equation Aipp = —Airpp. 

The Oseen's like term in Eq.Q provides a model de- 
scription of the hydrogen bonds network deformation en- 
ergy in the long wavelength limit. The value of the phe- 
nomenological parameter C « 10~^^cra^ , as well as the 
specific form and characteristic values of the function 
V (s^) describing the elastic energy of water polarization 
should be extracted from the experimental properties of 
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Figure 1: Typical polarization of the water molecules between 
lipid membranes shown in green. Both the charges of the lipid 
polar groups and the polarization of the water molecules are 
schematically represented by the corresponding charge sym- 
bols. 

the liquid as described in in |23]. Typically the liquid 
polarization is small, s ^ 1, and "the equation of state" 
function takes the usual Ginzburg-Landau form 

V{s')^p + Bs\ (2) 

where A = An/ie-l) « 0.16, e « 80 is the (large) 
dielectric constant, and B ~ 1. The parameter A 
characterizes the long-range interactions of the molec- 
ular dipoles, depends strongly on the temperature: A oc 
{T-Tc)/Tc 122, and vanishes &t T = Tc ~ 22%K 
roughly at the A— transition point in supercooled water 
|16l [TTl [181 [19] . On the contrary, the parameter B de- 
pends on the short-range physics and thus is practically 
independent of the temperature. The last two terms in 
Eq. ([ij describe the long-range dipole-dipole interaction 
of molecules, and interaction energy of the liquid in the 
external electric field, produced by charges that reside 
inside or outside the liquid. 

The cell membrane consists primarily of a thin layer 
of amphipathic phospholipids which spontaneously ar- 
range so that the hydrophobic "tail" regions are shielded 
from the surrounding polar fiuid, causing the more hy- 
drophilic "head" regions to associate with the cytoso- 
lic and extracellular faces of the resulting bilayer, as 
shown on Fig. ([iJ. Typically the width of the lipid 
bilayer is up to ^nm. Consider two membranes of the 
cross section area S separated by the water filled layer 
of width h from each other. Typically the membrane 
surfaces are comprised by molecular groups with no net 
charge, though with a considerable dipole moment. The 
lipid "heads" are hydrophilic, therefore water molecular 
dipoles arrange themselves in the direction of the lipid 
wall: s(±/i/2) — ±so, so ~ 1 (of course, < sq < 1: the 
value of So depends on the dipole moments of lipid polar 
groups). 

Since there is no external electric field in our system, 
Ep = — 47rP = — 47rPoS and the expression for the free 
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energy of the liquid reads 



C /ds\- 



G = +Uis)j, (3) 

where s = is the only non- vanishing component of the 
average water molecules polarization and U{s) = V{s'^) + 
2tts\ 

Variational minimization of the free energy (jlj with 
respect to the variations of the function s(r), SG = 0, 
gives the equation 



CcP^ 



+ {V'{3^) + 27r) s = 



(4) 



similar to the equation of motion of a particle of the mass 
C at the position s moving in one-dimensional quasi- 
potential U{s). Accordingly, the derivative ds/dz plays 
the role of the particle velocity velocity and the variable 
z serves as the time. The solution is a well known Id 
soliton 



dz 



± 



C 



where the exact value of the "energy" is a constant to 
be found by matching the boundary conditions at the 
membrane surfaces s(±/i/2) = sq: 



h 



/2C 



ds 



yr+iTM' 



(5) 



The functione(/i) simplifies in the large h limit: the main 
contribution to the integral comes the small s region 
where the integral diverges logarithmically, the function 
V{s^) can be neglected altogether, U{s^) « 27rs^, and 
therefore 



27rsQ exp(— /i/A), 



(6) 



where A = <^ h is the characteristic size of the 

soliton, A ^ Rjj w 0. 15-^0. 25nTO, the characteristic water 
orientation domain size, first introduced in |22l |23] . 

Eq. (§ can be used to transform the energy of the 
liquid layer to 



G 



2C 



h/2 



dz 



— eh 



/2C / " ds^/e + U{s^) - eh. 
Jo 



According to the standard definition the pressure, P 
— S~^dG/dh, can be expressed as: 



P 



'2CP. 



. de 



So 



ds 



de. 



The integral in the r.h.s. can be evaluated using Eq. ([5) 
and the asymptotic expression ([gJ for e(/i), so that the 
pressure is given by 

P = Piexp{-h/X), 

where the prefactor Pi = 27rPgSQ. This is exactly the 
dependence observed in experiments [^. The pressure 
is indeed positive, the membranes do repel each other, 
and the forces vanish exponentially quickly as soon as 
the membranes depart further apart than ft, ~ A. 

The physical reason of repulsion is the formation of 
polarization charges in the middle of the water layer due 
to inhomogeneous water molecules alignment. The po- 
larization charges are all of the same sign (positive for 
the example presented on FigjlJ and thus generate the 
electrostatic repulsion, which together with the hydro- 
gen bonds network deformation energy is ultimately re- 
sponsible for the repulsion of the membranes. For the 
maximum possible surface polarization, sq = 1, which 
corresponds to the large dipole moments of the lipid po- 
lar groups limit, we can estimate the maximum value of 
the prefactor: (-Pi)„£i^ ~ 3 • 10^° dyne/cm?. The value 
of the pressure obtained in this way agrees well with the 
results of the measurements reported for different types 
of the membranes: 5 • 10^ < {Pi)f.xp < 5 • 10"'^° dyne/cm^ 

m. 

Few concluding remarks should be added here. First, 
the repulsion pressure depends essentially only on sq = 
{Sz)-p, ^-component of unit vector S directed along the 
dipole moment of water molecule on the membrane inter- 
face r. This quantity decreases slightly as the tempera- 
ture increases in agreement with the molecular dynamics 
studies of |8]. This observation lets us put some more 
weight in favor of the polarization |9l UHl E] over the so 
called "entropy" hypothesis [12] of the repulsion pressure. 
Indeed, the polarization pressure P decreases slowly as 
the temperature rises, whereas any entropy-related effect 
should lead to the sharp increase of P, as explained in 

Second, the model defined by Eq. ([ll is very similar 
to the non-linear screening model introduced in |26l [27] 
and was originally applied for the calculation of the elec- 
trostatic forces in water. The non-linear screening model 
does not contain the Oseen term responsible for the short- 
range ordering of the water molecules. The scale A is the 
characteristic size of the water molecules domain (clus- 
ter) depends on C and thus can only appear in the com- 
plete model ([ij. This makes our model apparently the 
minimal continuous model capable of predicting repul- 
sion of hydrophilic membranes. We note that the sep- 
arate problem interactions of hydrophobic objects has 
been also extensively studied within another class of two- 
scale continuous water models p4, 31, 32 . We leave the 
research on possible convergence of the approaches and 
the relation between the scales of the models for a future 
publication. 

The work was supported by Quantum Pharmaceuti- 
cals. Phenomenological water models are used in Quan- 
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turn software to compute solvation energies, protein-drug vironments. 
interactions and drug transport properties in aqueous en- 
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